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1. Introduction 

Growing network models with preferential attachment, such as the Barabasi- Albert 
(BA) model [1], have become increasingly popular. They bring a natural explanation for 
the scalefreeness observed in complex networks, either natural or man made [2, 3, 4, 5, 6]: 
most networks exhibit a broad degree distribution falling off as a power law {fk ~ k"^). 
These model systems also provide a natural setting to investigate other, more refined, 
features of networks and other random growing structures. In this work we focus our 
attention onto leaders and lead changes. Luczak and Erdos [7] already describe as 'a kind 
of a race' the growth of connected components in the Erdos-Renyi model for random 
graphs [8], and refer to the largest component as the 'leader'. This terminology was 
then introduced in the physics literature by Krapivsky and Redner [9], in concomitance 
with the growing interest in social networks. As suggested by these authors, since the 
degree of a node may quantify the wealth of a corporation or the popularity of a person, 
it is natural to investigate questions such as: 

How does the identity of the leader (the node with largest degree) change in the course 
of time? What is the probability that a leader retains the lead as a function of time? 

The present work is devoted to an investigation of questions of this kind in the 
Bianconi-Barabasi (BB) fitness model [10, 11]. This model has a natural temperature 
parameter T. It interpolates between the BA model [1] (T = oo) and the record-driven 
(RD) growth model [12] (T = 0). It exhibits a condensation transition at some finite 
critical temperature T^, for an appropriate choice of the distribution of fitnesses. Our 
main goal is to extend to the full BB model the analysis of the statistics of leaders and 
lead changes, respectively performed in [12] and [9, 13] for the RD and BA models. 
We shall consider the case where the model has a condensation transition, in order 
to investigate the differences in behavior between the fiuid phase (T > T^) and the 
condensed phase (T < Tc). Besides the interest in leaders per se, the present study 
will also shed new light onto the dynamics of the condensed phase of the BB model, a 
subject so far only tackled in [14]. 

The BB model is a growing network model with preferential attachment. A new 
node enters the network at each time step. The node entering at time n is labeled by 
its index n. It attaches to a single earlier node {i — 1, . . . ,n — 1) with probability 

Z(n-l) • ^'-'^ 
The attachment probabilities Pn,i combine both a fit-get-richer feature, through the 
intrinsic fitnesses r]i of the nodes, and a rich-get-richer feature, through their dynamical 

degrees ki{n — 1). More precisely: 

(i) The fitness r]i is assigned to node i once for all and modeled as a quenched random 
variable. Fitnesses are conveniently parametrized as activated variables: 

7], = exp{-ei/T), (1.2) 

temperature T providing a measure of the ruggedness of the fitness landscape, whereas 
the activation energies Si are independent and identically distributed (i.i.d.) random 
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variables, drawn from some temperature-independent distribution. Nodes i with a low 
energy have a high fitness rji, and therefore a larger chance to attract new connections. 
The lower bound of energies is set for convenience to e = 0. The behavior of the density 
p{e) of the distribution for £ — )> turns out to play a central role in the model. We 
assume that p{6) starts as a power law: 

p{e)-e'-' {e^O), (1.3) 
with exponent ^ > 0. In numerical simulations we shall use the triangular distribution 

p{e) = 2£ (0 < £ < 1), (1.4) 
which has 9 — 2. 

(ii) The degree ki{n — 1) is the number of connections of node i just before node n 
enters the network. Nodes with a high degree, i.e., which are already well connected, 
have a larger chance to attract new connections. 

(iii) Finally, the partition function Z(n — 1) in the denominator of (1.1) ensures that 
the attachment probabilities are normalized: 

n 

Z{n) = J2mhin). (1.5) 

i=l 

As done in our previous studies [13, 15], we will hereafter consider two different 
initial conditions in parallel: 

Case A. Node 1 appears at time n = 1 with degree ki{l) = 0. At time n = 2 node 2 
attaches to node 1, so that ki{2) = k2{2) = 1. At time n = 3 node 3 can attach cither 
to node 1 or to node 2. Hereafter we make the choice of attaching node 3 to node 1, 
obtaining thus /ci(3) = 2, whereas ^2(3) = A;3(3) = 1. 

Case B. Node 1 appears at time n—1 with degree A;i(l) = 1. This amounts to saying 
that the first node is connected to a root, which does not belong to the network. It is 
natural to represent this connection by half a link. At time n — 2 node 2 attaches to 
node 1. We thus have ki{2) = 2 and k2{2) = 1. 

In both cases the network has the topology of a tree. The sum of the node degrees 
at time n equals twice the number of links L{n) in the network: 

f2Hn)^2Lin), (1.6) 

i=l 

with 2L^^\n) = 2n — 2 and 2L^^\n) = 2n — 1. Here and in the following, the 
superscripts (A) and (B) denote a result which holds for a prescribed initial condition, 
i.e.. Case A or Case B. We have therefore {k)^^\n) = 2 — 2/n and {k)^^\n) = 2 — 
which both yield 

(k) = 2 (1.7) 

in the thermodynamic limit, as expected for a tree topology. 

Owing to its temperature parameter T, the BB model realizes a continuous 
interpolation between two limiting models. 
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• The Barabasi- Albert (BA) model [1], for T = oo. 

In this limit we have r]i = 1 for all nodes i, so that node fitnesses do not play 
any role. The BA model is thus recovered. The attachment probability Pn,i is 
proportional to the degree ki{n — 1). The partition function reads Z{n) = 2L{n) 
(see (1.6)). 

• The record-driven (RD) growth model [12], for T = 0. 

In this limit, the growth of the network is record-driven. The current record node, 
i.e., the node whose fitness is the largest, attracts all the new connections, until it 
is outdone by the next record. 

Let us finally introduce some definitions used in the present work. 

(a) Leaders and lead changes 

A co-leader at time n is any node whose degree is equal to the largest degree 



The leader at time n is the node among the co-leaders whose degree reached the value 
kmiL^{n) first. We denote by I{n) the index of the leader at time n. Initial conditions are 
such that the first node is the first leader. There is a lead change at time n if the leader 
at time n is different from that at time n — 1. We call a lead any period of time between 
two lead changes. We denote by jO,{n) the number of leads up to time n. In other words, 
the number of lead changes up to time n is £(n) — 1. Some lead changes bring to the 
lead a node that has already been the leader in the past, whereas some other changes 
promote a newcomer. We denote by D{n) the number of distinct leaders up to time n. 
We also define the lead persistence probability S{n) as the probability that there is a 
single leader up to time n. 

(b) Records 

The record node at time n is the node whose energy (fitness) is the smallest (largest) 
met up to time n. A node is said to be a record if it belongs to the series of records, 

1. e., if it has been the record during some time lapse in the past. We denote by n(n) 
the probability that the leader at time n is a record. 

The setup of this paper is as follows. In Section 2 we revisit the continuum mean- 
field-like approach of [10, 11], emphasizing in particular finite-size effects near the critical 
point. Section 3 is devoted to the history of leaders and lead changes during a typical 
instance of the growth of the network, and Section 4 to an investigation of the degree 
statistics in the condensed phase. A brief discussion followed by a summary arc given 
in Section 5. An example of extreme-value statistics for i.n.i.d. random variables and 
a few complements on the record-driven growth process are respectively exposed in 
Appendix A and Appendix B. 

2. The continuum formalism revisited 

The main features of the model, and in particular the possible occurrence of a 
condensation transition, have been studied in the original works of Bianconi and 
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Barabasi [10, 11]. They use a continuum mean-field- like formalism, where both the 
discrete structure of the model and all fluctuations are neglected. In this section, we 
revisit this analysis in order to unveil some features of the model untouched so far, such 
as the importance of finite-size effects near the critical point, or the possible occurrence 
of an infinite hierarchy of condensates in the low temperature phase. 

2.1. Phase diagram 

Our starting point is the following exact recursion equation for the mean degree: 

{h{n)) - {h{n - 1)) = {pn,i) {n > i), (2.1) 

where brackets denote an average over the history of the network up to time n, i.e., 
both over the stochastic growth process and over the quenched random fitnesses. In the 
spirit of the continuum formalism, i.e., neglecting all fiuctuations and treating time n 
as continuous, we approximate the above equation as (see (1.1)): 

dkj, ^ i]iki 

dn ^ Z{n) ■ ^ ' ' 

We furthermore estimate the partition function in the denominator as Z[n) ~ Cn, with 

C = {rjk) (2.3) 

(see (1.5)). The differential equation (2.2) can then be integrated as 

h{n)^(^) . (2.4) 



This result is the cornerstone of the continuum approach. 

The mean degree ^(77) of a node with fitness rj can then be estimated by uniformly 
averaging (2.4) over the index i: 

where the rightmost expression is obtained by replacing the sum by an integral. The 
mean degree (k) can then be estimated by averaging (2.5) over the distribution of the 
fitness T). We thus obtain 

{k)^l + K{T,C), (2.6) 

and similarly 

{r]k) = CK{T,C), (2.7) 

with [10, 11] 

(see (1.2)). The equalities (1.7) and (2.3) consistently yield the condition 

K{T,C)^1, (2.9) 
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which determines the parameter C as a function of temperature. 

As said carUer, the main feature of the BB model is the possible occurrence of a 
condensation transition at some finite temperature Tc [10, 11], according to the value of 
the exponent 6 of the distribution of the energies £j (see (1-3)). 

• If 6' < 1, we have K(T, C) ^ oo as C — )■ 1"*" for any finite temperature. The 
condition (2.9) defines a regular function C{T), decreasing from C(oo) = 2 to C(0) = 1. 

• li e > 1, the integral K{T,C) is convergent as C 1+. The limit K{T,1) is a 
decreasing function of temperature, such that K{T, 1) — > oo as T — > oo and K{T, 1) — > 
as T — > 0. There is a unique critical temperature Tc such that 

i^(T„l) = l. (2.10) 

The condition ^ > 1 is somewhat similar to the condition d > 2 for the occurrence of a 
Bose-Einstein condensation in a free Bose gas. 

• In the high-temperature (fiuid) phase {T > Tc), the function C{T) is regular and 
decreasing from C(oo) = 2 to C{Tc) = 1 at the critical point. 

-k In the low-temperature (condensed) phase {T < Tc), C(T) = 1 stays equal to its 
critical value, whereas the r.h.s. of (2.6) keeps varying with temperature. This implies 
that a macroscopic fraction of the degrees escapes the sum rules (1.7) and (2.3). The 
analogy with the Bose-Einstein mechanism suggests that a condensation takes place [11]. 
Assume that some finite number J of nodes, labeled for j = 1, . . . , J, arc condensates, 
in the sense of having macroscopic degrees ki^ ~ R{j)n, proportional to the total number 
of nodes n, and fitnesses rji. !^ 1 close to optimal. Let F — J2j R{j) be the total size 
(i.e., degree fraction) of these condensates. Equations (2.6) and (2.7) read 

(k) = 1 + K{T,C) + F = 2, {r]k) = CK{T,C) + F = C, (2.11) 

hence C = 1 and 

F ^1- K{T,1). (2.12) 

The total size F of the condensates is thus predicted to be a well-defined function 
of temperature, increasing smoothly as temperature is lowered, from F = at the 
critical point to F = 1 in the T — )> limit. Let us emphasize that the continuum 
approach neither predicts the number J of condensates nor the individual sizes 
of the condensates. The possibility of there being more than one condensate in the 
low-temperature phase of the BB model is hinted at in [6]. We shall investigate these 
matters in Section 4. 

2.2. Finite-size effects 

We now turn to an investigation of finite-size effects near the critical temperature. Let 
us have a closer look at the sum appearing in (2.5). It is of the form 

n 

Sn{S) = j:t'-\ (2.13) 

i=l 
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where the exponent 5 = 1 — ?7/Cisin the range < (5 < 1. For fixed 5 > 0, we 
can estimate Sn{S) at large n by replacing the sum by an integral. We thus obtain 
Sn{S) = / 5 + ■ ■ ■, where the dots stand for a (5-dcpcndcnt correction of order unity. 
For 5 = 0, 5'„(0) is equal to the harmonic number of order n, and therefore diverges 
logarithmically as 5'n(0) Inn + 7b, where 7e is Euler's constant. In the scaling region 
where 5 is small and n is large, the sum behaves as 

p(51nn 1 

Sr,{5) « ^—^ . (2.14) 

This expression interpolates between the two limiting forms seen above. It also puts 
forward a time-dependent scale of 5, 

6. ^ r^, (2.15) 
mn 

to which we shall come back in Section 3. Inserting the scaling form (2.14) into (2.5), 
we obtain a finite-size correction to (2.6) of the form 

{k)^l + K{T, C) - AK{T, C), AK{T, C) = ' ^^'^^^ 

In the high-temperature phase (T > Tc), we have C{T) > 1, so that AK{T,C) ~ 
j^-t^(T) fg^iig Qg g^g pQwer law, with a positive correction exponent 

w(r) = i-i/c(r). (2.17) 

Finite-size corrections for generic quantities obey the same power-law fall-off. The well- 
known n~^^^ scaling of finite-size corrections in the BA model is recovered in the T — oo 
limit, since C(oo) = 2. 

At the critical point, we have C{Tc) = 1, so that the correction exponent lj{Tc) 
vanishes, and a more accurate treatment is needed. Assuming C = 1 exactly, and hence 
T Tc, we have l — rj/C ~ e/Tc. As a consequence, the expression (2.16) for AK(Tc, 1) 
is dominated by small values of the activation energy e, and simplifies to 

AK(T,,1)~AT, r£^-2e-(-i-")/^=d£ = ^^^^-;^^. (2.18) 
Jo (Inn)''"^ 

In order to investigate finite-size scaling around the critical point, we define an 
effective time-dependent critical temperature Tc^es{n) = Tc + AT(n) by requiring that 
C — 1, i.e., more precisely, that the mean partition function equals {Z{n)) — n at 
T — Tc^es{n). The shift AT'(n) provides a measure of the size of the critical region for a 
large but finite time n. It can be estimated by setting (k) = 2 in (2.16). We thus obtain 

^ , , AX(Te,l) 1 

The denominator 

dK(T,l) 1 ep(e)de 
is indeed convergent whenever 9 > 1. 
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Figure 1. Effective temperature Tc^cff('T-) against 1/ Inn. Symbols: numerical data for 
both initial conditions. Dashed line: extrapolation of the data for initial condition B. 



We now give a few explicit formulas in the case of the triangular distribution (1.4). 
This example of a distribution with 6 = 2 will be used in numerical simulations in the 
following. We have 



so that 



and 



KiT, C) 



KiT, 1) 



2e de 



2TE 



T - (m + T)e-'"/'^ 



m>l 



2TE 



(m + T)e-"*/^ 



m>l 



2tt^T 2 



T 



E 

»n>l 



2mT + 2T^)e 



-m/T 



(2.21) 



(2.22) 



(2.23) 



All the above series are exponentially convergent. The condition (2.10) yields 

Tc ^ 0.711 716 308, /i(Tc) ^ 1.896 474 189. (2.24) 

Figure 1 shows numerical results illustrating the slow logarithmic fall-off of the 
size of the critical region predicted by (2.19). Here and throughout the following, 
data are produced by means of a direct numerical simulation of the BB model for 
the triangular distribution (1.4). Data are averaged over 10^ independent realizations 
for each initial condition (A or B). The figure shows the effective temperature Tc cfr(^) 
against 1/ Inn. The observed linear convergence corroborates the prediction (2.19). The 
fit of the data for Case B (dashed line) extrapolates to Tc ~ 0.711, thus reproducing the 
critical point (2.24) to three digits. 
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3. The history of leaders and lead changes 

3.1. The index of the leader: a variational approach 

In order to estimate the index I{n) of the leader at time n, we start from the 
prediction (2.4) of the continuum formahsm for the degree of node i at time n, that 
we rewrite as 



Within this framework, which neglects the noise involved in the growth process, 
estimating the degree /Cmax(?^) and the index I{n) of the leader amounts to studying the 
largest value Zmax of the n random numbers Zi defined above and the index, hereafter 
denoted by /inid, where this largest value is reached. 

In other words, we are facing a problem of extreme-value statistics for independent 
but not identically distributed (i.n.i.d.) random variables. The classical extreme- 
value theory for independent and identically distributed (i.i.d.) variables is well 
documented [16]. In contrast, little is known for i.n.i.d. random variables (see [17, 
Section 4.2] for a recent review). The present problem provides an example of strongly 
i.n.i.d. variables. The support of the random variables, i.e., < < ln(n/i), indeed 
shrinks rapidly as a function of the index i. The index /inid is therefore expected to be 
localized in a relatively small range w <^n of values of i. 

In order to estimate the localization range w, we propose a variational approach, 
designed for the study of strongly i.n.i.d. random variables, i.e., whose supports have 
a strong dependence on their index. The predictions of the variational approach will 
be tested against exact results for a specific example in Appendix A. This approach 
proceeds as follows. Neglecting all numerical prefactors, we evaluate the index /inid 
as w. The minimal energy Smin reached in this range is evaluated as the smallest in 
a set of w i.i.d. variables Sj. According to a well-known argument of extreme-value 
statistics, e„un is such that the probability of having e < £min, i-e., £min (^^^ (1-3)), is of 
order 1/w. Putting both above estimates together, we obtain 



This expression has a non-trivial maximum as a function of w. If w is too large, the 
typical weight ln{n/w) is too small, whereas if w is too small, e^i^ is too large because 
of an insufficient sampling. It is natural to estimate the localization range w as the 
location of this maximum. We thus obtain the following implicit equation for w: 



Consistently neglecting prefactors and subleading corrections, we are left with the 
following prediction, both for the index Jinid where the largest value ^max is reached, 
and for the index /(n) of the leader in the original problem: 





(3.1) 




(3.2) 




(3.3) 




(3.4) 
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Figure 2. Black and red: mean index {I{n)) of the leader against Inn for T — 0.4 and 
both initial conditions. Blue: mean index (/inid) corresponding to the variables Zi. 
Black dashed curve: localization range w predicted by (3.3), translated horizontally 
for readability. 



In turn .Z'^ax and fc^ax can be estimated as 




These scaling predictions are the main results of this section. Let us note that the 
corresponding minimal energy scales as emin ~ T / Inn. The associated value of the gap 
5 = 1 — r^max ~ ^min/T, i.e., 5 ~ 1/ lun, is consistent with the scale (2.15). 

Before analyzing the consequences of the above results, in Figure 2 we present a 
check of the variational approach against numerical data for a temperature T = 0.4 in 
the condensed phase. The black and red curves show data for the mean index {I{n)) 
of the leader against ln?7, for initial conditions A and B. Both essentially yield the same 
results, except at short times. The blue curve shows data for the mean index (/mid) 
corresponding to the largest of the i.n.i.d. variables Zi introduced in (3.1). The data 
have been obtained by averaging over 10^ independent samples of these variables for each 
value of n. The black dashed line shows the localization range w as predicted by the 
implicit equation (3.3). For clarity, n has been multiplied by an arbitrary factor (i.e.. Inn 
has been translated). The good agreement between {I in)) and (/mid) demonstrates that 
the study of the variables Zi dictated by the continuum approach captures the essential 
features of the growth of the degrees. The qualitative agreement of the localization 
range w with all the data corroborates the validity of the variational approach. 

The main prediction of the variational approach is that the typical index of the 
leader at time n grows as I{n) ~ (Inn)^ (see (3.4)). This slow growth is expected to 
hold irrespectively of temperature T, i.e., both in the fluid phase (T > Tc) and in the 
condensed phase (T < Tc). 
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Figure 3. Black and red: ratio n/(fcmax) against Inn at the critical point {T — Tc), 
for both initial conditions. Blue: fit by a quadratic polynomial, slightly translated for 
readability. 



At variance with the index of the leader, the growth law (3.5) of the degree fcmax(^) 
of the leader depends on the phase of the model through the dependence of C itself on 
temperature. 

• In the fluid phase (T > Tc), we obtain the subextensive power-law growth 
^max ~ n^^'^^'^\ up to a logarithmic correction. The exponent 1/C(T) = 1 — a;(T) 
coincides with the maximal growth exponent of the expression (2.4). It is a 
decreasing function of temperature, varying continuously between the BA exponent 
1 /2 in the T — oo limit and the limit value 1 as T — )■ . 

• Right at the critical temperature (T = T^), the prediction (3.5) reads 

~ (3.6) 

The degree of the leader is therefore only logarithmically subextensive. The validity 
of this prediction is demonstrated in Figure 3, showing the ratio n/{kraax) against 
Inn. The data for both initial conditions are well represented by a quadratic 
polynomial (blue line), obtained by fitting the data for initial condition B, and 
slightly translated for readability. 

• In the condensed phase (T < Tc), the expression (3.5) is of little use, because of 
the presence of large finite-size corrections to C(T). The degree statistics in the 
condensed phase will be investigated in Section 4. 

3.2. A typical history of leaders and lead changes 



The scenario we are proposing for the typical history of leaders and lead changes is 
based on the following two ingredients. 
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Figure 4. Asymptotic probability I{oo{T) for tlie leader to belong to the series of 
records, against T/{T + Tc). Dashed lines: extrapolation of fits to the data. Symbols 
on right axis: corresponding quantities in the BA model with both initial conditions. 

First and most importantly, the index of the leader at time n grows as I{n) ~ (In n)^ 
(see (3.4)), irrespectively of temperature. There is therefore an asymptotic decoupling 
of time scales between n and I{n). 

The second ingredient is the high probability that the leader is a record. The 
probability n(n, T) that the leader at time n is a record, in the sense that it belongs 
to the series of record nodes, is observed to converge to an asymptotic value Yl^{T), 
irrespectively of the initial condition. This limit value is plotted in Figure 4 against 
T/(T + Tc). The crossovers to the BA and RD models at high and low temperatures 
(see the next two sections) prevent one from having accurate data at either end of the 
temperature range. The asymptotic probability is observed to decrease as a function 
of temperature, from IIoo — )■ 1 in the T — )■ limit, through ^.^{Tc) ~ 0.866, to a limit 
IIoo ~ 0.80 as T — 7- oo. The symbols to the right show the corresponding quantities 
for the BA model, i.e., ^ 0.613 and «i 0.703. The discrepancy between these 
numbers and the previous limit is a manifestation of the discontinuous behavior to be 
described in Section 3.3. 

Summarizing what we have learned so far about the leader in a large network 
consisting of n nodes, the index of the leader typically grows as (Inra)^, whereas it has 
a high chance of being a record. Therefore, as a rule of thumb: The Leader of Today is 
the Record of Old. The scenario for the history of leaders and lead changes to which we 
come up is essentially insensitive to the value of temperature, and is therefore expected 
to hold both in the fluid phase (T > Tc) and in the condensed phase (T < Tc). 

Several consequences can be drawn from this scenario. The typical number D{n) 
of distinct leaders up to time n is expected to grow (at most) as the typical number of 




■ BA (A) 

■ BA(B) 
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records up to time I{n) ~ (Inn)^, i.e., D{n) ~ ln7(n), or 

~ ^Inlnn. (3.7) 

The data for -D(n), to be presented in Figure 5, are indeed observed to grow very slowly 
with time n. A quantitative check of the asymptotic doubly logarithmic dependence 
is however practically impossible. The lead persistence probability S{n), i.e., the 
probability that the first node has kept the lead up to time n, is then expected to 
fall off exponentially with D{ri), i.e., as 

S{n) ~ l/(lnn)^ (3.8) 

3.3. High-temperature regime: crossover to the Barahdsi- Albert (BA) model 

In the T — > oo limit, the BB model becomes the BA model. For any finite time n, every 
quantity in the BB model goes continuously to the analogous quantity in the BA model. 
Some asymptotic long-time quantities may however exhibit discontinuous behavior. In 
technical words, the n — )> oo and T — )■ oo limits need not commute. An example is the 
asymptotic probability Hoc (7") for the leader to be a record, shown in Figure 4. The 
crossover scale can be determined by considering the index of the leader. In the BA 
model, the mean index of the leader has finite asymptotic values which depend on the 
initial state [13], namely (Z)^-' ~ 3.40 and (Z)^'' ~ 2.67. In the BB model, the prefactor 
of the scahng law (3.4) vanishes in the T — > oo hmit. The crossover in the n-T plane 
takes place for a time n ~ tba(Z^), with 

lnrBA(T)~T (T ^ oo). (3.9) 

For a fixed high temperature T ^ 1, the BB model behaves essentially as the BA model 
for n ^ T-BP^{T), whereas the asymptotic behavior of the BB model only appears for 
n > tba{T). 

3.4. Low-temperature regime: crossover to the record-driven (RD) growth process 

In the T — > limit, the BB model becomes the record-driven (RD) growth process [12]. 
In the latter model, the current record node attracts all the new connections, until 
it is outdone by the next record. Here again, some asymptotic long-time quantities 
may exhibit discontinuous behavior. The crossover scale can again be estimated by 
considering the index of the leader. The following fine of reasoning holds irrespectively 
of the presence {9 > 1) or of the absence {9 < 1) of a condensation transition at a 
finite critical temperature T^. In the RD growth process, as a consequence of the scale 
invariance of the underlying record process, the mean index of the leader grows linearly 
with time n, as {I{n)) ^ ^n. The analytic expression of the amplitude x ~ 0.275 765 is 
derived in Appendix B (see (B.ll)), using the formalism of [12]. In the BB model, the 
prefactor of the scaling law (3.4) diverges as T — )■ 0. The crossover time trb{T) is such 
that {{lnT-R£,)/Ty ~ trd, hence, to leading order, 

rKB{T)^T-' (T^O). (3.10) 




Figure 5. Data for the BB model for three low values of temperature T and initial 
condition B (colors) and for the RD process (black). Left: mean index {I{n)) of the 
leader against n. Right: mean number {D{n)) of distinct leaders against Inn. The 
dashed lines have the theoretical asymptotic slopes x and w. 



For a fixed low temperature T ^ 1, the BB model behaves essentially as the RD process 
for n <^ T-K£){T), whereas the true asymptotic behavior only appears for n ^ T-^-oiT). 

Figure 5 gives an illustration of the crossover of the BB model to the RD model 
on two characteristic quantities. The left panel shows {I{n)) against n. The RD data 
closely follow the linear law with slope x (dashed line), whereas the BB data bend 
toward the logarithmic growth (3.4). The right panel shows {D{n)) against Inn. The 
RD data closely follow the linear law with slope u ~ 0.624 330 (dashed line) (see (B.8)), 
whereas the BB data bend toward the doubly logarithmic growth (3.7). 

From an operational viewpoint, for a fixed temperature T, the crossover time can 
be defined as being the time for which (/(n, T)) = (/(n, 0))/2 is half its value in the 
RD growth process. Figure 6 shows times the crossover time tiit){T) thus defined, 
against temperature T, for both initial conditions. The measured intercept yields the 
behavior rRD(r) ^ 0.25/T^, in agreement with (3.10). 

4. Degree statistics in the condensed phase: an infinite hierarchy of 
fluctuating condensates 

The goal of this section is to give a picture for the dynamics of the condensates in the 
low-temperature phase (T < Tc). To do so, we investigate the statistics of the ordered 
sequence of largest degrees. For a given large time n, we rank the n nodes of the network 
in order of decreasing degrees: 

fc(l) > ^2) > /C(3) > • ■ ■ , (4.1) 

and we set 

kij) = Rij)^- (4-2) 




0.05 0.1 0.15 0.2 

T 



Figure 6. Plot of times ttie crossover time trd(2^) against temperature T for both 
initial conditions. Dashed line: extrapolation of a common fit to the data, yielding the 
intercept 0.25. 

The reduced variable Rq) represents the 'size' or 'weight' of the j-th largest degree, i.e., 
the degree fraction carried by the corresponding node. In particular A;(i) = fcmax is the 
degree of the leader at time n, and is the degree fraction carried by this leader. 

Let us recall that, according to the continuum approach of [10, 11], the total size of 
all the condensates is equal to the following deterministic (i.e., non- fluctuating) quantity: 

^ = E%) = 1-^(^'1)- (4-3) 

j 

(see (2.12)). As already discussed in Section 2.1, the continuum approach does not say 
anything about the number of condensates, nor about the mean values or the fluctuations 
of the individual sizes of the condensates. 

A first qualitative picture is given in Figure 7, showing five typical tracks of the 
degree fraction = k^^x/n of the leader against Inn. The left panel shows data 
in the RD growth process (T = 0), whereas the right one shows data at temperature 
T = 0.2. In both situations the degree of the leader exhibits large fluctuations around 
its mean value, that is (-R(i)) = w ~ 0.624 330 (see (B.8)) for the RD growth process, 
and (-R{i)) ~ 0.46 for T = 0.2. In the RD growth process, because of the scale invariance 
of the underlying record process, we observe a stationary signal on a logarithmic scale. 
At finite temperature, the tracks first resemble those of the RD process, however with 
some noise on a short time scale. The dynamics then slows down drastically when time 
exceeds a few times the crossover time trd (see (3.10)). We have Iutrd ~ 3.41 for 
T = 0.2 and initial condition A, whereas the rather sudden slowing down is observed 
near Inn ^ 4.5. Let us however note that fluctuations in the degree of the leader do 
persist forever, on a very slow scale, though. The number of maxima or minima is 
indeed expected to grow proportionally to the number of distinct leaders, i.e., as In Inn 
(see (3.7)). 
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Figure 7. Five typical tracks = fcmax/"- against Inn for initial condition A. Left: 
RD growth process (i.e., T = limit). Right: T = 0.2. 
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Figure 8. Mean values and standard deviations of the sizes of the larger two 
'condensates' (j = 1,2), and mean value (F), at temperature T = 0.2, against 
l/(lnn)^. Thin (resp. thick) lines show data for initial condition A (resp. B). 



Turning to a more quantitative viewpoint, we present in Figure 8 the mean 
values (-R(j)) and standard deviations a{R(^j)) = {{Rfj)) — (-R(i))^)^^^ of the larger two 
condensates (j = 1,2), again at temperature T = 0.2. The mean value {Y) of the 
moment variable Y, to be defined below (see (4.4)), is also plotted. The data have 
the following limiting values, irrespectively of the initial condition: (-R(i)) ~ 0.461, 
a(i?(i)) ^ 0.213, (i?(2)) ^ 0.100, a(i?(2)) ^ 0.075, and (Y) ^ 0.280. The choice l/(lnn)2 
of the abscissa is motivated theoretically by the prediction (3.6) for the largest degree at 
the critical point. A phenomenological analysis of the data, e.g. by fitting a power-law 
fall-off, would essentially yield the same asymptotic values. 
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Several conclusions can be drawn from the above results. The extrapolated values 
of (-R(i)) and {R{2)) are clearly different from zero, indicating that the condensed phase 
of the BB model is characterized by the presence of more than one condensate. The 
widths are comparable to the means {R(j)), demonstrating that the sizes 

are non-self-averaging. 

The above observations point toward the scenario of an infinite hierarchically 
ordered sequence of fiuctuating condensates. This situation is thus qualitatively similar 
to that of the RD growth process describing the model at zero temperature, recalled in 
Appendix B. Namely, the reduced variables it!(j) (j = 1,2,...) have a non-trivial joint 
distribution in the limit of a large network (n — )■ oo). In other words, the condensed 
phase is characterized by the presence of an infinite hierarchy of condensates, whose 
sizes R(j) are non-self-averaging but keep fiuctuating forever. The mean sizes of the 
successive condensates however fall off rapidly with the rank index j. 

One convenient tool to investigate this collection of fiuctuating condensates consists 
in introducing the second-moment variable 

which converges to the asymptotic variable 

r = E4)- (4-5) 

i>i 

The second-moment variable Y (and its higher-order generalizations Y^) have been 
introduced seemingly for the first time by Derrida and Flyvbjerg [18], in an investigation 
of the sequence of weights obtained by the iterative random breaking of an interval. The 
data plotted in Figure 8 (blue fines) suggest the existence of such a fiuctuating limiting 
variable Y, such that (Y) ^ 0.280 at T = 0.2. 

The dependence of the distribution of the hierarchy of condensates on temperature 
is illustrated in Figure 9. The theoretically predicted total condensed fraction F 
(see (2.12)) is plotted against reduced temperature T/T^ (thick black line), as well 
as data for the mean degree fraction of the leader {R{i)) (red symbols) and for the 
mean moment variable {Y) (blue symbols). All over the condensed phase, the mean 
degree fraction (-R(i)) of the leader is observed to be significantly smaller than the 
total condensed fraction. This demonstrates that the multipficity of condensates is a 
significant phenomenon over the whole condensed phase. 

We now turn to a more detailed study of the two endpoints of the condensed phase, 
corresponding respectively to T — )■ and T T^. The low-temperature (T 0) regime 
exhibits a crossover to the RD growth process, as already underlined in Section 3.4. The 
large symbols on the vertical axis of Figure 9 show the corresponding values in the RD 
growth process, namely (R) = u ^ 0.624 330 (see (B.8)) and (Y) = 1/2 (see (B.13)). 
These numbers are clearly far below the extrapolated T — >■ limits {R(i)) ~ 0.865 and 
(Y) a; 0.82, thus providing explicit examples of the discontinuous behavior announced 
in Section 3.4. The manifestation of the RD to BB crossover on the degree of the leader 




On leaders and condensates in a growing network 18 

1 

0.8 
0.6 
0.4 
0.2 



0.2 0.4 0.6 0.8 1 

T/T^ 

Figure 9. Black line: theoretically predicted condensed fraction F against reduced 
temperature T/T^. Symbols: numerical data for the mean degree fraction of the 
leader (red) and the mean moment variable (Y) (blue). Symbols on left axis: 

corresponding quantities in the RD process. 
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Figure 10. Mean degree fraction of the leader (i?(i)) — {k^a^) /n against Inn, for the 
BB model with initial condition A at several low temperatures (colors) and for the RD 
process (black). Dashed line: limiting value uj for the RD process. 

is illustrated in Figure 10. For a fixed low temperature T, the mean degree fraction of 
the leader, (-R(i)) = {kmax) /n, first follows the behavior characteristic of the RD process, 
up to a time of order the crossover time rRD(T) (see (3.10)), tending to converge to the 
RD limiting value u (dashed line), before it bends away from u and eventually saturates 
to a much higher value, close to the T — )■ limit {R(i)) ~ 0.865. 

In the vicinity of the critical temperature (T T^), the total condensed fraction 
F ^ yu(Tc)AT vanishes linearly with AT = — T. The corresponding slope /i(Tc) has 
been evaluated in (2.24). The mean degree fraction of the leader is observed to 





On leaders and condensates in a growing network 



19 



fall off much faster than linearly as AT — )■ 0, thus indicating that the leader carries 
a smaller and smaller part of the total size of the condensates as the critical point is 
approached. This observation can be corroborated by means of the following line of 
reasoning. Right at the critical point (T = T^), the degree of the leader is predicted 
to grow as (3.6). The corresponding degree fraction therefore vanishes logarithmically 
slowly, as ~ l/(lnn)^. Now, with the natural assumption that the degree of the 
leader obeys a finite-size scaling law around the critical point, the above time dependence 
at T = Tc can be turned, using (2.19), to the estimate 

~ (AT)^/(^-i). (4.6) 

For 9 = 2 we obtain the quadratic behavior ~ (AT)^, in agreement with the data 
shown in Figure 9. 

The scaling law (4.6) can be rephrased in the following form. The effective number 
of condensates, defined in an operational way as the ratio J — F/ {R{i)), diverges as 

J ~ (AT)-V(«-i) (4.7) 

in the vicinity of the critical temperature (T — )■ T~). The exponents of the power 
laws (4.6) and (4.7) diverge as 6* ^ 1, where the condensation transition disappears. 

5. Discussion 

After [12, 13, 15], the present work closes up a cycle of four papers devoted to the 
degree statistics in growing networks with preferential attachment, the main emphasis 
being on leaders and lead changes. We have focussed most of our attention onto the 
Bianconi-Barabasi (BB) fitness model [10, 11] and on its infinite-temperature and zero- 
temperature limits, respectively corresponding to the Bianconi-Barabasi (BA) model [1] 
and the record-driven (RD) growth model [12]. Very few works had been devoted to 
leaders in growing structures so far [7, 9, 19, 20]. 

Our main findings concerning the history of leaders can be summed up as follows. 

(i) At infinite temperature (BA model), the leader statistics keeps a memory of the 
initial state of the network. A typical history of the network involves finitely many 
distinct leaders, chosen among the oldest nodes [13]. The mean index of the leader, the 
mean number of distinct leaders, and the lead survival probability have finite asymptotic 
values, which depend on the initial condition. 

(ii) At zero temperature (RD process), the history of leaders inherits the temporal 
self-similarity of the underlying record process [12]. This property manifests itself as 
stationarity on a logarithmic scale. There are typically In n records and uj \nn distinct 
leaders up to time n. The degree and the index of the leader grow linearly with time, as 
^max ~ and / ~ Sn, where R and S keep fiuctuating, and admit non-trivial limiting 
distributions so that (R) = u ^ 0.624 330, and {S) = x ^ 0.275 765. 

(iii) In the generic finite-temperature situation, the main features of the history of 
leaders hold irrespectively of temperature, i.e., both in the fluid phase (T > T^) and 
in the condensed phase (T < Tc). The leader is a record node with high probability. 
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Quantity 


BA {T = oo) 


BB (T finite) 


RD (T = 0) 


{D{n)) 


{D)oo finite 


9 In Inn 


a; Inn 


inn)) 


(7)00 finite 


{{lnn)/Ty 




(A;max(n)) 




T>T,: nVC(^) 
T^T^: n/(\nny 


uin 



Table 1. Asymptotic temporal behavior of the key quantities: mean number {D{n)) 
of distinct leaders, mean index {I{n)) of the leader, mean degree (A:max(")) of the 
leader, in the various regimes of the model. Only the last observable is sensitive to the 
presence of a condensation transition at a flnite temperature T^. 

whereas its index grows on a logarithmic scale, as I ~ ((In n)/T)^. The number of 
distinct leaders up to time n therefore grows forever, although extremely slowly, on a 
doubly logarithmic scale: D{n) ~ In Inn. This picture can be captured into the motto: 
The Leader of Today is the Record of Old. 

The above scenario for the history of leaders relies on one essential ingredient, which 
was missing in the work by Fcrrctti and Bianconi [14], namely the separation of time 
scales induced by the very slow growth of the index of the leader. This explains why the 
conclusions of the present work are in many respects more precise than those of [14]. 

When temperature is either very high, or very low, the generic behavior of the BB 
model is eventually observed beyond a crossover time scale tba{T) which diverges at 
high temperature (whereas BA behavior is observed for shorter times) or trd(7') which 
diverges at low temperature (whereas RD behavior is observed for shorter times). As 
a consequence, histories ending with a single final leader never occur except in the BA 
model, i.e., strictly at infinite temperature. 

The above discussion can be recast in the language of the rcnormalization 
group. The BA and RD models respectively appear as infinite-temperature and zero- 
temperature unstable fixed points, whereas the BB model corresponds to a stable 
finite-temperature fixed point. The associated renormalization-group fiow is sketched in 
Figure 11. The situation is the exact opposite of that usually met in critical phenomena. 
Furthermore, the fixed point describing the critical behavior around the condensation 
transition temperature T^ is not visible in this flow. Table 1 summarizes the asymptotic 
behavior of three key quantities (number of distinct leaders, index of the leader, degree 
of the leader) in the various regimes of the model. 
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All the predictions recalled so far arc insensitive to the presence of a condensation 
transition at a finite temperature T^. Only the last observable (/Cmax(?^)) depends on 
the position of temperature T with respect to Tc. If the BB model has no condensation 
transition (^^ < 1), the subextensive growth law (/Cmax(^)) ~ n^^^^'^^ holds for all positive 
temperatures. 

Another motivation for the present work was the desire of a better understanding 
of the dynamics in the condensed phase {T < Tc). On the one hand, the analogy 
with the Bose- Einstein mechanism put forward originally in [11] suggests a single static 
condensate with degree k^a:^{n) ~ Fn, whose relative size F is predicted by the mean- 
field-like approach (see (2.12)), and therefore does not fiuctuate in time. On the other 
hand, for some other stochastic processes leading to condensation, such as the zero- 
range process, the condensate is known to perform some ergodic stationary-state motion 
over the system, on a time scale which diverges as a power of the system size (spatial 
extent) [21]. The situation of growing structures is made more complex by the fact that 
the system size is time itself, so that strictly speaking there is no stationary state. 

At variance with earlier studies [11, 14], the present work leads to the picture of an 
infinite hierarchy of condensates, whose degrees grow linearly with time, as k^) ~ R{j)n. 
The relative sizes it!(j), generahzing the reduced degree of the leader — R — kmax/n, 
sum up to the predicted static condensed fraction F. The individual sizes of 
the condensates are however non-self-averaging. They rather keep fluctuating forever, 
and admit non-trivial temperature-dependent limiting distributions over the whole low- 
temperature condensed phase. 
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Appendix A. Extreme-value statistics for i.n.i.d. variables: an example 

In this Appendix we investigate the sequence of i.n.i.d. random variables: 

Zi^iUi (i = l,...,n), (A.l) 

where the Ui are i.i.d. variables uniformly distributed on the interval [0,1]. We are 
mostly interested in the distribution of the largest value .^max among the n random 
numbers Zi, and of the index Jinid where this largest value is reached. Our goal is to 
test the variational approach proposed in Section 3.1, by comparing its predictions to 
the exact results derived below. 

A.l. Variational approach 

In the present situation, the variational approach proceeds as follows. The random 
variable Zi lies in the interval [0,i], whose upper bound i takes its largest value for 
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i = n. The index /inid is therefore expected to be locahzed in a relatively small range 
of width w around n, i.e., n — /inid ~ w. On the other hand, the largest value t/max 
among the w variables Ui sampled in that range is such that the probability of having 
U > f/max) i-e., 1 — ?7max) IS of Order 1/w. Neglecting every numerical prefactor, the 
estimates n — /inid ~ w and 1 — ?7max ~ 1/w yield 

Tl 

n-Z^^r^w + —. (A.2) 
w 

This expression has a non-trivial minimum for w — y/n, where both terms are equal. 
The variational approach thus predicts the following scaling relations for large n: 

n - ^max ~ 2(n - /inid) ~ V^- (A. 3) 

A.2. Exact results for finite n 

The example (A.l) is simple enough to allow one to derive exact results on the 
distribution of ^max and /inid, for a finite number n of variables. 

Let us start with the largest value ^max- The common distribution function of the 
uniform variables Ui on [0, 1] reads 

f for < 0, 
u for < 1* < 1, (A.4) 
1 for u> 1. 



Fu{u) = Prob{C/i <u}^ < 



The distribution function of Zmax is therefore 



Fz^M = Prob{Zn,^ < ^} = n Prob{^[/i <z] = Y[Fu[-). (A.5) 

1=1 i=i ^ ^ ^ 

Setting Int z — n — k, the latter expression reads explicitly 

i=n-k+l * 

The mean value of ^max is given by 

{Zra^)= / F',_{z)zdz= / (1-Fz_(^))dz. (A.7) 

v 

Using the expression (A. 6), we obtain 

(^max)=n-5:^-^ / zUz 

" (n ~ A;)! {n - k + - {n - kf+^ 

_ " (n-fc)! (n-fc)^+^ 

n! (A; + l)(A; + 2)- ^ " ^ 

Let us now turn to the index /inid- Consider a fixed index i = 1, . . . ,n and a fixed 
value Zi of the random variable Z^. For any j ^ i, the probability to have Zj < Zi, i.e.. 
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n 


1 


2 


3 


4 


5 


6 


7 


(•^max) 


1 
2 


13 
12 


125 
'72 


3503 
1440 


5687 
1800 


49253 
12600 


19797847 
4233600 


(-^inid) 


1 


7 


31 


55 


1033 


829 


61153 


4 


12 


16 


240 


160 


10080 



Table Al. Mean largest value (.^max) and mean index (Jinid), as given by the exact 
results (A.8) and (A. 12), for n up to 7. 



Uj < Zi/j, reads Fu{zi/ j). The conditional probability that Zi is the largest therefore 
reads 



Prob{Z, = Z^^\i, Zi}^l[Fu[-]. (A.9) 



Setting Int Zi = n — k, wc have exphcitly 

Prob{Z, = Z^^\i, Zi} = ilLzil! i (A.IO) 

n\ 

The distribution Pj = Prob{/inid = index Jinid can be obtained by integrating 

the latter expression against the distribution of Zf. 

1 n 

Pi^ - PT0h{Zi = ZmaxK Zi} dZi 

I Jo 

n (^_^)! .n-k+1 

= E ^ — / 

k=n-i+l •^"-'^ 

" {n-k)\ (n - A; + 1)^ - (n - A;)^ 

k=n—i+l 

The mean value of the index li^id is given by 

n 

(/inid) = E ^ 
i=l 

-It^-^^iin-k + ir-in-kr) {2n + l-k). 

^ k=\ 

-^E^^(--e (A.12) 



2 .t; n\ 



Table Al gives the first few values of (^max) and (/inid), as given by the exact rational 
expressions (A.8) and (A.12). In spite of the simplicity of these numbers, which calls 
for a combinatorial interpretation, we have not been able to relate them to anything 
known in the literature. 
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A. 3. Asymptotic results for large n 

Let us now turn to the analysis of the large n behavior of the above results. The 
expression (A. 6) for the distribution function of Zmax can be estimated as 

nl 

where we have used the approximation z ~ Int z — n — k, as well as Stirling's formula. 
The above result shows that n — Z^^x. scales as A; ~ ^pn. More precisely, setting 

^max = n - Xy/n, (A.14) 
and differentiating (A.13), we obtain the following distribution for the scahng variable X: 

fx{x) = xe-""'/^. (A. 15) 

Along the same lines, the expression (A. 11) for the distribution of /inid can be recast as 

P,^- r e-'^'/^^^MA;. (A.16) 

n J n—i 

This result shows that n — /iuid again scales as A; ~ \/ri. Setting 

/inid ^n- Y^, (A. 17) 

the expression (A.16) yields the distribution of Y: 

fy{y) = /f e-^'^'do; = /I erfc (A.18) 



The mean values of the scaling variables read {X) — ^Jnji and (Y) — ^7r/8, so that 

{n - Z^ax) ~ 2(n - /inid) ^ sj^. (A. 19) 

This result fully corroborates the outcome (A. 3) of the variational approach. The 
predicted factor 2 exactly holds for the mean values, as we have {X) — 2(Y). 
More generally, higher moments of the scaling variables are related through {X'^) — 

(p+l)(n- 

Finally, corrections to the scaling formulas for the distributions of ^max and /inid 
can also be extracted by a more careful analysis of (A. 8) and (A. 12). Skipping details, 
we just give the expansions for the mean values, including subleading constant terms: 
, V /Trn 7 , -r ^ /Trn 5 ,. 

(^max) ~ ^ - V T ^ 6' ^ '"'^'^ ~ ~ V T ^ 6' ^ 

Figure Al shows the differences n — {Z^^ and n — (/inid) against y^. The straight 
blue hues show the asymptotic expansions (A.20). 




Figure Al. Plot of the differences n — {Z„^ax) and n — (/inid) against ^/n. Straight 
blue lines: asymptotic expansions (A. 20). 



Appendix B. Some complements on the record-driven growth process 

In this Appendix we give a few useful complements on the record-driven (RD) growth 
process [12]. In the RD process, only the nodes which belong to the series of records grow 
by attracting new connections. We label the record nodes by their rank m = 1,2,... 
in the series of records. We use the asymptotically exact description of the theory of 
records in terms of continuous variables. The m-th record node is born at time Nm 
(m = 1, 2, . . .), with Al = 1, and with the recursion 

= (B.l) 

where the Um are uniform i.i.d. random variables between and 1. 

Throughout the following it will be sufficient to consider the stroboscoped process 
at record times n = Nm- The distribution of all dimensionless quantities of interest can 
indeed be shown [12] to be the same at a record time A^ and at a generic time n. 

B.l. Distribution of the degree and of the index of the leader 

The leader at time n = N^n is one of the earlier record nodes, i.e., node A^ for some 
/i = l,...,m — 1. In the language of [12], /x is a leader, whereas the subsequent ones 
(/i + 1, . . . , m — 1, if any) are non- leaders. We have therefore 

n = Nm, I{n) = A^, = A^m+i - ^m- (B-2) 

The key point of the analysis consists in introducing the reduced variables 

p _ fcmax(r^) ^ _ /(n) 

rim — ; — ) \^-'^) 

n n 
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and in observing that they obey the recursion relations 
< Um+i < ^ 



1 + Rm 
^ < Urr,,+^ < 1 




(B.4) 



l + Rr, 

These equations can be viewed as defining a stochastic dynamical system with co- 
ordinates Rm and Sm subjected to a time-dependent noise Um- The dynamical system 
has a unique invariant measure. It therefore makes sense to consider the hmiting 
variables R and and especially their mean values {R) — uj and (S) = x, so that 

(/^max(^)) ~ wn, (/(^)) ^ xn- (B.5) 

The variables Rm are autonomous, in the sense that the recursion for the Rm only 
involves themselves and the noise Um- The variables Rm are asymptotically distributed 
according to the distribution of a variable R with density Jr between and 1. The 
Laplace transform of the density reads [12]: 

/i?(5) = (e-^/^) = l-e-^(^), (B.6) 

with 

E{s) = f (B.7) 
In particular the mean value of R reads (R) = u, with 

POD 

/ e-'-^('Ms fa 0.624 329 988. (B.8) 

^0 

This number is known as the Golomb-Dickman constant. It appeared in the framework 
of the decomposition of an integer into its prime factors [22], and in the study of the 
longest cycle in a random permutation [23, 24, 25]. 

At variance with the Rm, the variables Sm are not autonomous. The recursion (B.4) 
for the Sm indeed involves, besides themselves and the noise Um, the variables Rm as 
well, through the conditioning inequalities. The study of the variables Sm is therefore 
more intricate. In the body of this paper we are mostly interested in the mean value 
X — (S) . This quantity can be evaluated by means of the formahsm of [12] . Hereafter 
we only give the main hues of the derivation, skipping details. We translate labels of the 
record nodes by /i, so that record node number is leader, whereas the subsequent ones, 
labeled k — 1, . . are non-leaders. This is precisely the framework of [12, Section 5.2]. 
The recursion (B.4) yields 

Ri^l- U,, S, = U,, 

i?2 = (1 - Ui)U2, S2 = UiU2, (B.9) 

i?3 = (1 - U,)U2U^, ^3 = U.U^Us, 

and so on, so that 

S, = ^^^=1^- (B.IO) 




Figure Bl. Densities fa and fs of the distributions of the variables R and S. 



The joint distribution of Ri and Rk is explicitly known in the form of a double Laplace 
transform [12, Eqs. (5.43), (5.44)]. Some algebra leads to the requested analytical 
expression for the mean value x = (S)'- 

X = ^J^ (1 -e-*)e-^(')(e-* + E(s) -^(2s))ds ^ 0.275 765 063. (B.ll) 

To close up, we present in Figure Bl the distributions fji and fs, respectively giving 
the asymptotic distribution of the ratios k^a,x{n)/n and I{n)/n for the RD growth 
process at a generic time n. These very accurate data are obtained by iterating the 
recursion relations (B.4) 10^° times. 

B.2. Distribution of the ordered sequence of large degrees 

We now turn to the distribution of the whole ordered sequence of large degrees, using 
the notations of Section 4. The known relationship between the statistics of records and 
that of cycles in a random permutation [23, 24, 25, 26, 27] tells us, among other things, 
that all the admit limiting distributions, generalizing the distribution of i? = 
recalled above and plotted in Figure Bl. 

Closed-form expressions for all the moments of all the reduced variables have 
been derived by Shepp and Lloyd [25]: 

The mean values (-R(j)) and standard deviations a{R(^j)) = ((i?^^-)) — (i?(j))^)^/^ of the first 
few variables are listed in Table Bl. These numbers demonstrate the strong hierarchical 
ordering of the degrees of the larger nodes. The mean values {R{j)) indeed fall off 
as . The number of visible large nodes therefore grows logarithmically with n. We 
thus recover the logarithmic law of the typical number of records up to time n. 
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j 


(^0)) 




1 


0.624 330 


0.192 144 


2 


0.209 581 


0.112 044 


3 


0.088 316 


0.067037 


4 


0.040 342 


0.039 797 



Table Bl. Mean value and standard deviation of the first few reduced variables -R(j), 
evaluated from the expressions (B.12) derived in [25]. 



The expression (B.12) implies that the mean value of the moment variable Y defined 
in (4.5) reads 

(y) = I- (B.13) 

This result can be recovered alternatively as follows. The values of the moment 
variable at record times have been shown [12, Eq. (6.6)] to obey the recursion 

Ym+l = Ul^^Ym + (1 - Um+l)\ (B.14) 

with the same notations as in (B.4). Taking the mean value of both sides, we 
recover (B.13). We obtain similarly (F^) = 7/24, {Y^) = 139/720, and so on. 
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